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A future large- volume liquid scintillator detector would provide a high-statistics measurement of 
terrestrial antineutrinos originating from /^-decays of the uranium and thorium chains. In addition, 
the forward displacement of the neutron in the detection reaction i> e + p — > n + e + provides 
directional information. We investigate the requirements on such detectors to distinguish between 
QJ , certain geophysical models on the basis of the angular dependence of the geoneutrino flux. Our 

analysis is based on a Monte-Carlo simulation with different levels of light yield, considering both 



^ ■ unloaded and gadolinium-loaded scintillators. We find that a 50 kt detector such as the proposed 

LENA (Low Energy Neutrino Astronomy) will detect deviations from isotropy of the geoneutrino 
flux significantly. However, with an unloaded scintillator the time needed for a useful discrimination 
between different geophysical models is too large if one uses the directional information alone. A 
Gd-loaded scintillator improves the situation considerably, although a 50 kt detector would still 
need several decades to distinguish between a geophysical reference model and one with a large 
neutrino source in the Earth's core. However, a high-statistics measurement of the total geoneutrino 
flux and its spectrum still provides an extremely useful glance at the Earth's interior. 



"Electronic address: |hochmuth@ph.t um.de| 
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I. INTRODUCTION 

The first measurement of the v e flux from natural radioactive elements in the Earth 
has triggered a lot of excitement about the future of "neutrino geology" 

1MB. 

The 

geoneutrino 1 flux could deliver new information about the interior of the Earth, in particular 
its radiochemical composition, and thus shed new light on the question of how the Earth 
and other planets formed. Such an ambitious programme requires detectors of the next 
generation that are able to provide much larger event rates. 

One possible future detector that may well serve this purpose is LENA (Low Energy 
Neutrino Astronomy) that has been proposed by several of us |(| for high-statistics solar 
neutrino spectroscopy, for spectroscopy of the cosmic diffuse supernova neutrino background 
(DSNB), as a detector for the next galactic supernova, and to search for proton decay Q. 
Present design studies for LENA assume 50 kt of liquid scintillator that would provide a 
geoneutrino rate of roughly one thousand events per year from the dominant 

v e + P — > n + e + (1) 

inverse beta-decay reaction. 

While liquid scintillator detectors do not provide direct angular information, indirectly 
one can retrieve directional information because the final-state neutron is displaced in the 
forward direction. The offset between the e + and the neutron-capture locations can be 
reconstructed, although with large uncertainties. The CHOOZ reactor neutrino experiment 
has demonstrated a nontrivial directional sensitivity in that it was able to locate the reactor 
within 18° half-cone aperture (68% C.L.) on the basis of 2500 reconstructed events in a Gd- 

n 

loaded scintillator |8J. Therefore, it is natural to study the requirements for a future large- 
volume liquid scintillator detector to discriminate between different geophysical models of the 
Earth that differ both by their total neutrino fluxes and the neutrino angular distributions. 

We will consider two types of detectors. Motivated by current design studies for LENA 
we will consider a 50 kt detector using a PXE based scintillator. However, in this case it is 
difficult to locate the neutron-capture event on protons because a single 2.2 MeV gamma is 
released that travels on average 22.4 cm before its first Compton interaction. Therefore, the 



1 It is understood that this term refers to electron-antineutrinos emitted by the Earth's natural radioactive 
elements. 
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event reconstruction is relatively poor. As a second case we study a Gd-loaded scintillator 
where the neutron capture releases on average 3 photons with a total of 8 MeV, allowing for 
a much better event reconstruction. Moreover, the spatial resolution can likely be improved 
beyond the original CHOOZ experiment if one uses photomultipliers that are fast enough 
(time jitter ~ 1 ns) to use the photon arrival times for the event reconstruction. 

For the geoneutrino flux we will consider a continental and an oceanic location. In each 
case we will use a reference model and exotic cases with an additional strong neutrino source 
in the Earth's core or with enhanced neutrino emission from the mantle. 

We will begin in Sec. [0] with a discussion of the principle of geoneutrino detection in 
large- volume scintillator detectors as well as possible scintillator properties. In Sec. II 111 we 
introduce our geophysical models. In Sec. HVI we turn to the main part of our work, a Monte- 
Carlo study of the discriminating power of different detectors between different geophysical 
models and conclude in Sec. El 



II. GEONEUTRINO DETECTION 



A. Directional information from neutron displacement 

In a scintillator detector, geoneutrinos are measured by the inverse beta-decay reaction 
Eq. (Q) with an energy threshold of 1.8 MeV. The cross section is 

a = 9.52 x lO" 44 cm 2 (2) 

MeV MeV v ; 

where E + is the total energy of the positron and p + its momentum. The visible energy 

E v - ls = E + + m e always exceeds 1 MeV because the positron annihilates with an electron 

of the target. By measuring the visible energy one can determine the neutrino energy as 

E v ~ E v - ls + 0.8 MeV because the kinetic energy of the neutron is typically around 10 keV 

and thus negligible. After thermalization the neutron is captured by a nucleus, thus tagging 

the inverse beta decay reaction. 

Kinematics implies that the neutron is scattered roughly in the forward direction with 

respect to the incoming neutrino this being the key ingredient for obtaining directional 

information. The maximum scattering angle is 



2E U A - (A 2 - ml) 



cos# max = 1 , (3) 



where A = m n — m p and m n , m p and m e are the masses of the neutron, proton and positron, 
respectively. In this extreme case, the neutron and positron momenta are perpendicular to 
each other. For the maximum relevant geoneutrino energy of 3.2 MeV one obtains cos 9 XDax = 
0.79, which is equivalent to # max = 37.8°. However, most geoneutrinos have energies much 
closer to threshold and the maximum angles are much closer to the forward direction. The 
average displacement between the neutron and positron events is then theoretically found 
to be about 1.7 cm j^|. 

The reactor experiment CHOOZ, using a Gd-loaded scintillator, has measured an average 
neutron displacement from the e + event of 1.9 ± 0.4 cm j^J. However, once the neutron has 
been thermalized by collisions with protons, it diffuses some distance before being captured 
so that the actual displacement varies by a large amount for individual events. In a PXE 
based scintillator the average time interval until capture on a proton is 180 /xs, leading to 
an uncertainty a of the displacement of about 4 cm for the x-, y- and z-direction P|. With 
Gd loading a is reduced to approximately 2.4 cm j^| because the neutron diffusion time is 
much shorter, on average about 30 /xs sj]. 

B. PXE-based scintillator 

One option for the proposed LENA detector is to use a scintillator based on PXE (phenyl- 
o-xylylethane, CieHig). PXE has a high light yield, it is non hazardous, has a relatively high 
flashpoint of 145°C, and a density of 0.985 g/cm 3 [l(3]. A possible admixture of dodecane 
(C12H26) increases the number of free protons and improves the optical properties. A blend 
of 20% PXE and 80% dodecane shows a decrease in light yield of about 20% relative to pure 
PXE, an attenuation length of about 11 m and an increase in the number of free protons 
by 25% Q. 

In this paper we consider a detector with a total volume of about 70 x 10 3 m 3 . This could 
be realized with a cylindrical detector of 100 m length and 30 m diameter. An outer water 
Cherenkov detector with a width of 2 m acts as a muon veto. In order to shield against 
external gamma and neutron radiation a fiducial volume of about 42 x 10 3 m 3 with a total 
number of 2.5 x 10 33 free protons as target can be realized using a scintillator mixture as 
mentioned above with 20% PXE and 80% dodecane. In Monte-Carlo calculations the light 
yield of events in LENA has been determined pj. For events in the central detector region 
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the yield N pe , measured in photo-electrons (pe) per MeV energy deposition, can be expressed 
as N pe ~ 400 pe/MeV x c, where c is the optical coverage which depends on the number and 
aperture of the photomultiplier tubes (PMTs). A maximal coverage c max ~ 0.75 can not be 
exceeded so that we assume the maximal light yield to be around 300 pe/MeV. For instance, 
the use of 12,000 PMTs with a diameter in aperture of 50 cm would result in an optical 
coverage of about 30% and a light yield N pe ~ 120 pe/MeV. This can be obtained either by 
using PMTs like in the Super-Kamiokande experiment or by smaller PMTs equipped with 
light concentrators as they were developed for the Counting Test Facility (CTF) at the Gran 
Sasso underground laboratory Q|. For events off the axis of the cylinder the light yield 
would be enhanced. Hence, low-energy spectroscopy even in the sub-MeV region should be 
possible in LENA. 

For a detection of the positron-neutron displacement the ability of the detector to locate 
the absorption position of both particles is crucial. The experimental reconstruction of 
both events is possible by analyzing the arrival times and the number of photons in each 
individual PMT. The position uncertainty depends on the total yield of registered photo- 
electrons. In the CTF, the measured position uncertainty was around 10 cm in each direction 
for events with 300 photo-electrons and it was shown that the uncertainty scales with the 
inverse square-root of that number [ljj. Therefore, we will use a Gaussian distribution for 
the uncertainty of the positron event reconstruction with equal width in each direction of 



where N pc is the light yield and E V - 1S the visible energy released by the positron. 

In PXE-based scintillators the neutron is captured by a proton with nearly 100% efficiency 
within an average time interval of about 180 /is, subsequently emitting a 2.2 MeV gamma. 
This photon has a mean free path of 22.4 cm before its first Compton scattering so that 
the event reconstruction is much more uncertain than for the positron event. We have 
simulated this case by taking into account multiple Compton scatterings of the 2.2 MeV 
gamma. The position of each gamma emission, representing the position of the neutron 
capture, is reconstructed by composing the energy-weighted sum of each Compton scattering 
event, taking into account the instrumental resolution. The distribution of the reconstructed 
position in each direction follows roughly a Lorentzian form. In Fig. ^ we show the radial 
distribution of the reconstructed positions of these events for light yields of N pe = 50, 300 and 
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FIG. 1: Monte-Carlo simulation of the radial distribution of a 2.2 MeV 7-quantum in an unloaded 
PXE scintillator. The curves are for light yields of N pe = 50, 300 and 700 pe/MeV as indicated. 

700 pe/MeV. Increasing the light yield does not significantly narrow the distribution because 
its width is dominated by the large Compton mean free path of the 2.2 MeV photon. With 
reduced light yield the position of the maximum as well as the mean value of the distribution 
shifts towards larger values. This is caused by the increased uncertainty of the instrumental 
resolution. 

C. Gadolinium-loaded scintillator 

In PXE-based scintillators the neutron is captured by a proton with nearly 100% efficiency 
within an average time interval of about 180 /xs, subsequently emitting a 2.2 MeV gamma. 
However, one can alsoput an additive into the scintillator to enhance the delayed-neutron 
signal. Both lithium Q and gadolinium 3 loaded scintillators have been used in this way 
for neutrino experiments or are planned in the future Q|. With a concentration of about 
0.1% by mass, Gd captures neutrons with 90% efficiency and isotropically emits a total 
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energy of about 8 MeV in a gamma cascade with an average of 3 photons. Note however, 
that in a Gd-loaded PXE/dodecane scintillator the light yield is reduced by typically 10-20% 
compared to an unloaded scintillator. 

For Gd-loaded scintillators we use the measurements performed in the CHOOZ exper- 
iment. The neutron response of the detector was measured with a 252 Cf source and the 
position uncertainty of a neutron capture event was 19 cm at la in each direction Q|. The 
distribution was Gaussian and the light yield was measured to be N pe = (125 ± 5) pe/MeV. 
The reconstruction was performed using only the information about the amplitude distribu- 
tion of the PMTs. Based on these numbers we assume the position resolution of the neutron 
event for the LENA detector to be Gaussian with a width in each direction of 



In CHOOZ no information of the arrival time was used for the position reconstruction. How- 
ever, the CTF measurements demonstrated that this yields the most valuable information, 
provided the time response of the PMTs is fast enough, i.e. a time jitter not much larger 
than 1 ns. Therefore, in a Gd-loaded scintillator one probably could achieve a much better 
neutron-event reconstruction so that our estimates are conservative. 

D. Backgrounds 

KamLAND has reported 152 events in the energy region relevant for geoneutrinos within 
a measuring time of 749 days and 3.5 x 10 31 target protons. From these events 127 ± 13 are 
dee to background Q. The m08 t relevant backed for the KamLAND site are reactor 
antineutrinos (80.4 ± 7.2 events). For the LENA detector positioned in the underground 
laboratory CUPP (Centre for Underground Physics in Pyhasalmi) in Finland (longitude: 
26° 2.709' E,latitude: 63° 39.579' N, 1450 m of rock (4060 m.w.e.)) this background would 
be reduced by a factor ~ 12, as the site is far away from reactors. Hence we expect for 
LENA at CUPP a reactor background rate of about 240 events per year in the relevant 
energy window from 1.8 MeV to 3.2 MeV. This background can be subtracted statistically 
using the information on the entire reactor neutrino spectrum up to ~ 8 MeV. 

Another important background for KamLAND is induced by radio impurities. A large 
concentration of the long-lived isotope 210 Pb is present in the KamLAND scintillator. In the 




(5) 
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decay chain of 210 Pb the a-emitting isotope 210 Po is present. Thus the reaction 13 C(a, n) 16 
can occur, mimicking the signature of geoneutrinos due to neutron scattering on protons and 
the subsequent neutron capture. The number of these background events in KamLAND is 
estimated to be 42 ± 11 [l7(. However, with an enhanced radiopurity of the scintillator, the 
background can be significantly reduced. Taking the radio purity levels of the CTF detector, 
where a 210 Po activity of 35±12/m 3 d in PXE has been observed Q], this background would 
be reduced by a factor of about 150 compared to KamLAND and would account to less than 
10 events per year in the LENA detector. 

An additional background that imitates the geoneutrino signal is due to 9 Li, which is 
produced by cosmic muons in spallation reactions with 12 C and decays in a /3-neutron cas- 
cade. Only a small part of the 9 Li decays falls into the energy window which is relevant 
for geoneutrinos. KamLAND estimates this background to be 0.30 ± 0.05 (l^. At CUPP 
the muon reaction rate would be reduced by a factor ~ 10 due to better shielding and this 
background rate should be at the negligible level of ~ 1 event per year in LENA. 



III. MODELS OF THE EARTH 



A detailed density profile of the Earth, the Preliminary Reference Earth Model, was con- 
structed by Dziewonski and Anderson in 1981 by monitoring seismic activities Based 
on the examination of meteorites and solar system materials (moon rocks and dust) and 
available Earth material, geologists have deduced a model for the distribution and concen- 
tration of elements in the Earth, the Bulk Silicate Earth Model [3]. The most dominant 
and abundant radioactive isotopes are 238 U, 232 Th and 40 K; their decays heat the Earth. The 
present-day total energy loss through the Earth's surface is about 40 TW or 82 mW/m 2 . 
The ratio of the energy production due to radioactivity to the total heat flow at the surface 
is known as the Urey ratio. In the Bulk Silicate Earth Model this ratio is assumed to be 
0.5, attributing 20 TW of the Earth's heat loss to radioactivity. Other estimates take the 
Urey ratio to be as large as 0.8 j3, 2ll |. 

The flux of geoneutrinos is directly linked to the rate of radioactive decays and to the 
generated heat. Therefore, it is of great interest to measure the geoneutrino flux and thus 
deduce the main contributor to the heat production. The elemental abundance ratios in the 
Bulk Silicate Earth Model are Th/U « 4 and K/U « 1.14 x 10 4 . According to this model, 
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radioactive isotopes are only in the crust and mantle because they are lithophile whereas the 
core is void of any significant amount of uranium, thorium or potassium. Different estimates 
for their abundances in the crust and mantle differ by factors of 2-3. An excellent overview 
of these abundances and their spread is given in the GERM Reservoir Data Base (2^]. 

An Italian group of physicists has, in cooperation with geologists, constructed a reference 
model for the abundance values of uranium, thorium and potassium. They used the values 
referenced in GERM and derived a mean value for each element 2^. We will implement 
these abundances into our Reference Model of the geoneutrino angular distribution. For 
a discussion of the angular spectra of this model and its uncertainties see Ref. Our 
Reference Model is in accordance with the Bulk Silicate Earth constraint of a heat production 
of 20 TW due to radioactive decays in the crust and mantle. 

Other authors have speculated about the presence of uranium and thorium in the core, for 



24j,|25j. Moreover, the Earth's magnetic field 



example in the context of a putative georeactor 
is not yet fully understood, but seems to be generated by a complex interaction between core 
and mantle. The core itself is too hot to sustain a permanent magnetic field. Therefore, it 
is assumed that the magnetic field is powered by a geodynamo, where in a simplified picture 
the magnetic field is generated by the motion of the liquid outer core. One problem of this 
picture is the unknown energy source for the geodynamo. Powering it by energy from the 
inner core leads to cooling and solidification. This process is constrained by the present-day 
size of the inner core, implying that it must be younger than 1.7 Gyr much less than the 
Earth's age of 4.5 Gyr. Therefore, several authors have concluded that there are radioactive 
elements in the core, providing heat and sustaining the geodynamo [26 1. 

The geoneutrino flux depends sensitively on location. The oceanic crust is depleted in 
radioactive elements, whereas the flux on the continents is dominated by the crust. Thus 
an experiment situated in the Pacific Ocean, e.g. on Hawaii, would have better access to the 
oceanic crust and the mantle. For an experiment located on a continent we have assumed 
a thickness of 50 km for the crust, implying a total neutrino flux in our Reference Model 
of 4.2 x 10 6 cm -2 sec -1 from uranium and 4.1 x 10 6 cm -2 sec" 1 from thorium decays. For 
an oceanic site we have chosen the crust to be rather thick (50 km), but not included any 
sediments. If one wanted to determine the mantle contribution, the oceanic crust would 
be a background to the measurement so that the assumption of a thick oceanic crust is 
conservative. The neutrino fluxes in this case are 1.25 x 10 6 cm -2 sec -1 from uranium and 
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TABLE I: Annual event rates for 2.5 x 10 33 target protons. Flavor oscillations have been included 

with a global reduction factor of 0.57. 

Model Continental Crust Oceanic Crust 

Reference Model 1.02 x 10 3 0.29xl0 3 

(A) 20 TW core 1.35xl0 3 0.62xl0 3 

(B) lOTWcore 1.19xl0 3 0.45xl0 3 

(C) 20 TW Lower Mantle 1.44xl0 3 0.71 xlO 3 

0.88 x 10 6 cm ~ 2 sec" 1 from thorium decays. The zenith-angle distributions of the reference 
geoneutrino flux corresponding to continental and oceanic locations are shown in Fig. |21 
Our assumption of a uniform oceanic crust with a thickness of 50 km has been made for 
computational reasons. Still, this model for the Hawaiian detector site is valid as our event 



rate is slightly lower than the rate found in 



In addition, as will become obvious in 



Sec. II VI changes in the crustal thickness of an oceanic site by even an order of magnitude 
can not be resolved by the detector. 

Besides our reference model we consider three "exotic" cases A, B and C, each of them 
either with a continental or an oceanic crust. 

(A) Fully radiogenic model with additional uranium and thorium in the core, accounting 
for 20 TW additional heat production. (Integrated neutrino flux increase of about 
32% relative to the reference model in a continental location, and 116% in an oceanic 
location.) 

(B) Same with 10 TW in the core. (Flux increase of 16% and 58%, respectively.) 

(C) 20 TW in the Lower Mantle. (Flux increase of 41% and 148%, respectively.) 

The zenith-angle distributions of the neutrino fluxes from these models have been determined 
along the lines of Ref. jj| and is shown in Fig. El 

To obtain the event rate in a scintillator detector, neutrino flavor oscillations have to be ac- 
counted for by including a global v e survival-probability factor of 0.57 as measured by Kam- 
LAND H- 

Matter effects for oscillations are not important because of the small geoneutrino 
energies. Moreover, for geoneutrino energies of 1.8-3.2 MeV and Am 2 = 7.9 x 10~ 5 eV 2 the 
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FIG. 2: Zenith-angle distribution of the geoneutrino flux where 8 = corresponds to the vertical 
direction. Shown is the total flux without flavor oscillations. Upper panel: Continental crust. Lower 
panel: Oceanic crust. 
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vacuum oscillation length is 57-101 km. Including distance and energy dependent survival 
probabilities is a negligible correction to a global reduction factor 27[. The annual event 
rates corresponding to our models, including the reduction factor, are shown in Tab. for a 
50 kton detector with a fiducial volume corresponding to 2.5 x 10 33 protons. 

Up to now we have assumed that the exotic heat source in the Earth's core is caused 
by uranium and thorium decays, i.e. the neutrino spectrum from this additional source was 
taken to be identical with the geoneutrino spectrum from the crust and mantle. However, 
the possibility of a natural reactor in the Earth's core ( "georeactor" ) has been discussed 

nn 

in the literature \24L 12 5| . In this case the neutrino flux could be similar to that from an 
ordinary power reactor with energies reaching up to about 8 MeV. With this assumption the 
total 4-7r-georeactor neutrino flux can be estimated to be ~ 1.9 x 10 23 s" 1 for a thermal 
power of 1 TW. Taking into account neutrino oscillations, the distance to the center of the 
Earth, and the detection cross section we calculate an event rate of about 210 y _1 TW -1 in 
LENA. At Pyhasalmi one would observe about 2200 events per year due to neutrinos from 
nuclear power plants. Assuming a systematic uncertainty for the neutrino flux from the 
power plants of 6.5%, as suggested in Q], we conclude that LENA will be able to identify a 
georeactor of > 2 TW after one year of measurement with a 3a significance. The influence of 
the uncertainty of the mixing angle 9u is negligible, because the flux of both the georeactor 
and the power plants depends on # 12 in the same way. 



IV. MONTE-CARLO STUDY 



To study the power of directional discrimination of a large liquid-scintillator detector 
we have performed a Monte-Carlo simulation of a large number of geoneutrino events and 
the corresponding directional reconstruction. We have assumed that the detector response 
is independent of the event location, i.e. only the spatial separation between the event 
v e +p n + e + and the location of neutron capture is relevant. However, as pointed out in 
chapter II, we consider a position resolution of point like events located at the central axis 
of the detector. As the light yield and hence also the position resolution increases for off- 
axis events our assumption is conservative. We have assumed that, on average, the neutron 
capture point is displaced by 1.9 cm in the forward direction relative to the e + event in 
agreement with the CHOOZ measurement Q]. Moreover, we have assumed that neutron 
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FIG. 3: Monte-Carlo example of reconstructed zenith-angle distribution of 25,000 events in a Gd- 
loaded detector. Upper panel: Neutrinos injected from the horizontal direction (cos# = 0). Lower 
panel: Neutrinos injected from the vertical direction (cos# = 1). 



diffusion before capture causes a Gaussian distribution around this mean value with a width 
a x = a y = <j z = 4.0 cm for an unloaded PXE-based scintillator, whereas 2.4 cm is taken for 
a Gd-loaded scintillator as described in Sec. HTJ 

In addition to this distribution, the main uncertainty originates from the reconstruction 
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of both events. For the positron event we have assumed that the reconstructed location 
follows a Gaussian distribution with a width given by Eq. (J3J). The actual spread of relevant 
visible energies is small so that we have always used E vis = 1.4 MeV as a typical value. For 
an unloaded scintillator, the reconstruction of the neutron event introduces an even larger 
uncertainty; we have used a distribution as in Fig. Q appropriate for the given light yield. 
For a Gd-loaded scintillator the uncertainty can be well fitted by a Gaussian distribution 
with a width given by Eq. 

For the sake of illustration we discuss the reconstructed zenith-angle distributions for 
the two extreme cases where all neutrinos come from the horizontal direction or all of them 
come vertically from below, i.e with cos# = 1. Figure El displays the reconstructed zenith- 
angle distribution of both cases, each generated with 25,000 neutrino events and using a 
Gd-loaded detector with a light yield N pe = 300 pe/MeV. We conclude that, given the 
relatively poor angular reconstruction capability of scintillator detectors, the only angular- 
distribution information that can be extracted is the slope of the distributions shown in 
Fig. El Put another way, one can extract the total event rate and the dipole contribution 
of the angular distribution, whereas a determination of higher multipoles is unrealistic. 
Therefore, we write the reconstructed zenith-angle distribution in the form 



configuration of geophysical model and detector type. 

The event rates for our fiducial detector size with 2.5 x 10 33 target protons and different 
geophysical models have already been reported in Tab.|U What remains to be determined by 
means of a Monte-Carlo simulation are the corresponding coefficients p and their uncertainty. 
In Tab. |H] we show the results for p for different cases, always assuming a light yield of 
300 pe/MeV. The uncertainty a v of the measured p value scales with the inverse square root 
of the number of events N so that s p = a p \fN is a quantity independent of N. The value 
of s p can be derived analytically for p — 0, yielding 



which is valid for all p <C 1. We have checked with our Monte Carlo that Eq. (J7J) indeed 
applies to all p values of interest to us. 




(6) 



where the event rate N and the coefficient p are the two numbers that characterize a given 




(7) 
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TABLE II: Coefficient p for the reconstructed zenith-angle distribution for different Earth models 
and different detector types, always assuming a light yield of 300 pe/MeV. 



Model 



Coefficient p for scintillator detectors 
Unloaded PXE Gd-loaded 



0.0377 
0.0521 
0.0459 
0.0485 

0.0646 
0.0824 
0.0772 



Continenal Crust 

Reference Model 0.0283 

(A) 20 TW core 0.0377 

(B) 10 TW core 0.0333 

(C) 20 TW Lower Mantle 0.0351 
Oceanic Crust 

Reference Model 0.0468 

(A) 20 TW core 0.0597 

(B) 10 TW core 0.0560 



The number of events it takes to distinguish at the la level between an isotropic event 
distribution (p = 0) and the actual coefficient is given by iV lCT = (s p /p) 2 = (3/4) p -2 . For 
our reference model at a continental site we find Ni a w 500 events, for an oceanic site about 
200 events. In order to distinguish a geophysical model % from model j at the la level, the 
required number of events is 

H 3 1 , (8) 

lfT (Pi-Pj) 2 2 (Pi-Pj) 2 ' 

A detection at the na level requires n 2 times more events. 

In the same way as for Tab. |H] we have calculated the slope p for different light yields 
of the scintillator and have determined the number of events it takes to distinguish each 
of the exotic models from the reference case. In Fig. 0] we display A r lcr for these cases and 
the continental- crust situation as a function of the light yield N pe , both for an unloaded 
PXE-type detector and a Gd-loaded one. In Fig. El we show the same for the oceanic crust. 
In the oceanic location we do not show model C because it corresponds to an increased flux 
from the mantle, i.e. it is essentially identical with the reference model except for the total 
flux and thus can not be distinguished on the basis of the zenith-angle distribution. 

Of course, the time required to achieve this discriminating power depends on the detector 
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size. For our fiducial volume with 2.5 x 10 33 target protons as in LENA one needs to 
scale with the event rates shown in Tab. HJ In a continental- crust location, all models 
produce an event rate of roughly 1000 events per year, in full agreement with the KamLAND 
measurement Assuming for example a realistic light yield of 120 pe/MeV for the 

LENA detector (unloaded) one would have to measure roughly 150 year to distinguish even 
our most optimistic model A from the reference model. Therefore, even with optimistic 
assumptions a 50 kt detector would need several decades for distinguishing in a meaningful 
way between different geophysical models on the basis of the angular event distribution. 
Moreover, detector backgrounds should be included in a realistic assessment. 



V. CONCLUSIONS 



A future large-volume scintillator detector such as the proposed 50 kt LENA would 
provide a high-statistics measurement of the geoneutrino flux. The event rate would depend 
strongly on the detector location, notably on whether an oceanic site such as Hawaii is 
chosen where a reference event rate of about 300 per year (50 kt scintillator) is expected 
or a continental site such as the Pyhasalmi mine in Finland where the reference rate would 
be about 1000 per year. Therefore, the geoneutrino flux could be measured with high 
significance and would allow one to distinguish between different geophysical models. 

The forward displacement of the neutron in the inverse beta decay detection reaction 
provides directional information on the geoneutrino flux. We have studied if this effect can 
be used to distinguish between different geophysical models, notably if one could diagnose 
a strong exotic energy source in the Earth's core under the assumption that its neutrino 
spectrum is identical with that emitted by the crust and mantle. While a deviation from 
an isotropic flux can be ascertained with high significance, we find that a 50 kt detector is 
too small to distinguish between different geophysical models on the basis of the directional 
information alone, except perhaps for extreme cases and optimistic assumptions about the 
detector performance. 

In our study we have only used the neutrino flux from the Earth, ignoring the contribution 
from power reactors because it depends strongly on location. For example, in Pyhasalmi 
the neutrino flux from power reactors adds roughly 25% to the counting rate in the energy 
window relevant for geoneutrinos. This contribution is not negligible, but it does not change 
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FIG. 4: Number of events needed to distinguish between models A, B or C and the continental- 
crust reference model at lcr significance. The points correspond to the values calculated with the 
Monte Carlo. Upper panel: Gd-loaded scintillator. Lower panel: Unloaded PXE-type scintillator. 
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FIG. 5: Same as Fig. 21 for oceanic crust. Model C was not included because the neutrino flux has 
almost the same angular distribution as the reference model so that almost infinitely many events 
are needed for a discrimination. 
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our overall conclusions. 

We have also estimated the sensitivity of a LENA type detector for determining a hy- 
pothetical georeactor in the Earth's core. As a possible location the CUPP underground 
laboratory in Pyhasalmi (Finland) was chosen and the background due to nuclear power 
plants was calculated. At CUPP a 2 TW georeactor could be identified at a statistical level 
of 3a after only one year of measurement. 

In summary, large- volume scintillator detectors of the next generation will be extremely 
useful to study the interior of the Earth in the "light of neutrinos". However, the prime 
information will be the total geoneutrino flux and its spectrum. It would be extremely chal- 
lenging to use the directional information alone to distinguish between different geophysical 
models. 
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